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The combined effect of frustration and correlation in electrons is a matter of considerable interest 
of late. In this context a Falicov-Kimball model on a triangular lattice with two localized states, 
relevant for certain correlated systems, is considered. Making use of the local symmetries of the 
model, our numerical study reveals a number of orbital ordered ground states, tuned by the small 
changes in parameters while quantum fluctuations between the localized and extended states produce 
homogeneous mixed valence. The inversion symmetry of the Hamiltonian is broken by most of 
these ordered states leading to orbitally driven ferroelectricity. We demonstrate that there is no 
spontaneous symmetry breaking when the ground state is inhomogeneous. The study could be 
relevant for frustrated systems like GdI-2, NaTiC>2 (in its low temperature C2/m phase) where two 
Mott localized states couple to a conduction band. 

PACS numbers: 



I. INTRODUCTION 

Geometric frustration in correlated systems brings 
about a variety of phenomena and is a major area of inter- 
est in the condensed matter community presently. Sys- 
tems having two dimensional (2D) layered structure with 
triangular lattice (for example, transition metal dichalco- 
genidesi~— , cobaltates^, Gdl^^and its doped variant^, 
NaTiO^ 1 ^, NaVOz^i etc.) are known to come up with 
a host of cooperative phenomena like valence and metal 
insulator transitions, charge, orbital and magnetic order, 
unconventional superconductivity, excitonic instability 1 
and possible non-Fermi liquid states&i£. These systems 
pose a challenge to theoretical understanding as the un- 
derlying geometric frustration of the triangular lattice, 
coupled with strong dynamic fluctuations, give rise to 
a large degeneracy at low temperatures and competing 
ground states close by in energy. A consequence of this is 
a fairly complex phase diagram^ and the presence of soft 
local modes strongly coupled to the itinerant electrons^. 

We motivate the model by looking at two systems Gdl 2 
and NaTi02 . Each layer of Gd ions in Gdl2 form a 2D tri- 
angular lattice. They are well separated from each other 
by intervening layers of large Iodine ions and do not in- 
teract significantly with each other. From recent band 
structure (LSDA) calculations^ it is known that three 
nearly degenerate, spin-polarized d-orbitals (d z 2, d x 2_ y 2 
and d xy ) cross the Fermi level. Calculations involving dy- 
namical local correlations^ show that these three nearly 
degenerate d-orbitals further break down to two doubly 
degenerate (d x 2_„2 and d X y) localized levels below and 
one extended (d z 2 ) level across the Fermi level. From ex- 
perimental studies 8 , it is observed that the ground state 
is insulating and is likely to be orbitally ordered in a three 
sublatticc fashion at half filling (each orbital occupying 
one sublattice) while a small doping away from half-filling 



leads to phase segregation. An effective Falicov-Kimball 
model 13 of spinless Fermions with two (degenerate) local- 
ized bands and one itinerant band is proposed recently 
for this system^. It would, therefore, be interesting to 
look for such orbitally ordered states in this context. 
The corresponding Hamiltonian, then, is 

(ij) £ja=l,2 

+u £ (/l/iadk)+^/E(^/ii4/i 2 )- (i) 

Here d\,di are, respectively, the creation and annihila- 
tion operators for electrons in the itinerant band and 
ftia > fia ar e the same for the two localized bands at the 
site i. The first term in Eq.(l) is the kinetic energy of 
d-electrons on a triangular lattice (only nearest-neighbor 
hopping is considered) while the second term represents 
the degenerate energy levels Ef of the /i and /2 elec- 
trons. The third term is the on-site Coulomb repulsion 
between d- and /-electrons. The last one is the local 
repulsion between the /-electrons. 

NaTi02 is another system where the above model can 
be applied. It also has a layered structure with alternat- 
ing NaO and TiO slabs and shows a correlation-driven 
metal- insulator transition 10 . The Ti 3+ ions with one 
electron in the d-orbitals are arranged on a triangular 
lattice and are octahedrally coordinated to oxygen ions. 
In the low temperature monoclinic phase the TiOg octa- 
hedron gets distorted, leading to a splitting in the triply 
degenerate t2 9 orbitals of Ti 3d: two orbitals are pushed 
below the Fermi energy while one straddles it. Orbital 
order due to electronic correlations has been predicted 
in this system earlier—. While the absence of magnetic 
order renders a consideration of spin degeneracy redun- 
dant, correlations in the 3c? band would further localize 
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FIG. 1: (colour online) /-electron configurations are shown at each site for (a) n/ 1 = n/ 2 = 1/8, (b) n/ 1 =n/ 2 = 1/4 and 
(c) nf t = rif 2 = 1/3. Black and green circles correspond to sites occupied by fi and /2-electrons respectively and open circles 
correspond to sites with no /-electron occupancy. 
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FIG. 2: (colour online) d-electron densities for different U 
and Uf values (set of four figures (a) to (d)) for nf ± = rif 2 = 

1/3. 



the two bands below Fermi level and the effective model 
above could as well describe the low energy dynamics of 
this system. 



The Hamiltonian conserves local /-electron occupation 
numbers n/,; Q = fj a fi a owing to the local U(l) gauge in- 
variance in the absence of f-d hybridization. Therefore, 
[nf t i a ,H] = and uii a — fj a fi a are good quantum num- 
bers taking values only or 1. The local conservation 
also implies that the Hamiltonian may be written as, 



H 



5> 

(ij) 



(ui)dfdj + Ef ^ + u f^2 WilWi2 ( 2 ) 



where h^(uj) = -ty + (U J2 a ^i.a ~ fi)Sij. 

We set the scale of energy as the nearest neighbor hop- 
ping t = 1. The eigenvalue spectrum of this Hamil- 
tonian, is easily obtained by numerical diagonalization 
on a finite size triangular lattice with periodic bound- 
ary condition. In order to calculate the average values 
of physical quantities, the classical Monte Carlo method 
using Metropolis algorithm can then be employed by 'an- 
nealing' over a subset of configurations of the "classical" 
variables {oJi. a }. This approach is known to give reliable 
results for the FKM, even on a triangular lattice with 
macroscopic degeneracies. The details of the method is 
reported elsewhere^. 

For the systems discussed above, the number of elec- 
trons in the d-band is one at every site in the undoped 
limit. Therefore, in the following calculations, we restrict 
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FIG. 3: (colour online) Number of sites doubly occupied by 
localized f-electrons versus U for different [// (versus Uf for 
different U (inset)) for nj 1 = n/ 2 = 1/8. 



ourselves to the half-filled limit, i.e. Nf + Nd — N where 
Nf, Nd and N are the total number of /-electrons, d- 
electrons and the number of sites in the lattice respec- 
tively (we have used N = 144 in all the calculations 
and checked a few results with larger N). In addition, 
Nfi — Nf2 — Nf/2, as the two /-levels are degenerate. 



II. RESULTS AND DISCUSSION 

In Fig.l we show the ground state configurations for 
different U and Uf at various / electron concentrations. 
Each column shows the progression of states from dis- 
order to order to disorder and finally order again as a 
function of U and Uf for different nf values. Fig. 1(a) 

presents the case for n/j (= ^p) = n/ 2 = 1/8. In the 
absence of any interaction between the two localized /- 
electron states (Uf — 0), no particular order is observed 
even on increasing the value of U. In fact, a finite U 
favors double occupancies of /-electrons, leaving a larger 
part of the lattice for (^-electron motion. Therefore, with 
increasing U there is an increasing tendency of phase 
separation in real space between localized and itinerant 
electrons. Keeping U small (for example, U = 1) and 
making Uf finite begins to remove the local double occu- 
pancies of / electrons, reduces the overall kinetic energy 
of d-electrons and leads to an ordered stripe-like pattern 
of the two localized electrons (/i and f%)- This stripe 
pattern appears at a higher Uf when U is raised. A tran- 
sition from a state of disordered doubly occupied sites to 
a stripe-like ordered pattern as a function of Uf at fixed 
U is shown in Fig.l. This is a discontinuous transition 
appearing at a critical value of Uf (which increases with 
U) and the ordered phase remains stable up to high val- 
ues of Uf. Note that on raising U / , the double occupancy 
is removed very quickly (see later) and the ground state 



is fairly insensitive to a rise in Uf further. Similar trends 
are seen for other fillings too (e.g., nf 1 — rif 2 = 1/4 and 
n h = n li = V^, shown in Fig. 1(b) and 1(c) respec- 
tively) , albeit with different real space patterns of the / 
electrons: for nf 1 = rif 2 = 1/4 it is a bi-stripe orbital or- 
der whereas for nf t = nf 2 = 1/3 it is a three sub-lattice 
structure involving /i , fi and empty (no /-electron) sites 
decorating the vertices of each triangle. 

One of the most interesting consequences of the var- 
ious orbital orders seen in Fig.l is that many of them 
break the inversion symmetry of the Hamiltonian. This 
will induce spontaneous displacements of the ions and 
the resulting lattice distortion will lead to ferroelectric- 
ity with a spontaneous finite polarisation. This is due to 
a specific realization of the unconventional orbital order 
that breaks inversion symmetry. It is driven by electronic 
correlation, very different from the conventional route to 
ferroelectricity. 

In Fig. 2 we present the d-electron densities at each 
site. It is interesting to observe how the competition 
among different interaction energies (such as kinetic en- 
ergy of the d-electrons, on-site correlation energies be- 
tween d- and f-electrons as well as between two localized 
f-electrons) affect the transport of the system. The d- 
electron density plots at each site (presented in Fig. 2 for 
7if t = nf 2 = 1/3) clearly show the appearance of metallic 
and insulating and re-entrant metallic phases by tuning 
the parameters U and Uf. From Figs. 2; (b) and (c) we 
see that on increasing U to a high value, keeping U / fixed, 
the system goes from a stripe like ordered configuration 
to a phase separated state. This transition is accompa- 
nied as well by an insulator to metal transition as we 
discuss below. Again by keeping U fixed and increas- 
ing U f , the ordered structure reappears and the system 
becomes insulating (Figs. 2; (c), (d)). By tuning the in- 
teraction parameters U and Uf one can switch between 
these metallic and insulating states. Note that in some 
cases there are also sites substantially occupied by elec- 
trons of all three types (/i, f% and d). Number of sites 
doubly occupied by localized f-electrons predominantly 
plays a crucialrole in determining the critical U for the 
metal-insulator transitions discussed above. We show the 
dependence of this double occupancy (do) on the inter- 
action parameters U and Uf for n/i = n/2 = 1/8 in 
Fig. 3. The double occupancy do is extremely sensitive to 
Uf, dropping off to zero rapidly for Uf as small as 0.25 
for various value of U and nf j(x . It increases with U as 
expected and saturates quite rapidly thereafter. These 
features are quite general and are obsereved at other fill- 
ings of n/ Q also. 

In order to investigate the electronic properties of these 
ordered phases further we have calculated the gapi£ in 
the spectrum for each of these phases. Fig. 4(a) shows 
the gap as a function of U / for a series of U values, for 
the filling 1/8. The disorder to order transition is indeed 
accompanied by a metal to insulator transition (MIT), 
the disordered phase has no gap at the the Fermi en- 
ergy, whereas the ordered states are insulating with a 
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FIG. 4: (colour online) Energy gap A (a) as a function of Uj for different values of U (b) as a function of U for different Uf at 
n h = »V 2 = 1/8- 



finite gap. In addition, the gap appears at a critical [// 
and remains constant beyond (independent of Uf as long 
as U is fixed). The discontinuous nature of the transi- 
tion is clearly visible in Fig. 4. On increasing U the same 
feature is observed, though the critical U / for MIT in- 
creases. As Uf is the correlation energy between two 
localized /-electrons and U is the same between local- 
ized and itinerant electrons, there is a competition be- 
tween these two interactions and therefore the critical 
Uf depends strongly on U. In Figs. 4(a) we show the 
energy gap as a function of Ut at fixed values of U for 
ri/i = Ti/2 = 1/8. An interesting observation can be 
gleaned from the energy gaps if we keep U f fixed (say 
Uf = 1) and vary U. For U = there is no gap in the 
spectrum and on increasing U, the gap increases almost 
linearly with U (Fig. 4(b)). Above a certain value of U 
(for Uf = 1, at around U = 3.2), however, the gap sud- 
denly drops to zero and the system becomes metallic. 
Therefore, as a function of U and at a fixed U / , there 
is a metallic phase appearing at large U. This seem- 
ingly counter-intuitive result is understood from the fact 
that there are two competing processes at hand (com- 
ing from U and Uf). At U = the d-electrons are free 
and the system is metallic. As U increases they avoid 
sites occupied by /-electrons and a gap appears. There 
is also a tendency towards f-d phase segregation as dis- 
cussed above. A finite value of Uf, however, strongly 
disfavors double occupancy (as in Fig. 3) and spreads the 
/-electrons out. Therefore, when, at a higher U, double 
occupancy (with /i and f% electrons) is energetically fa- 
vorable over a joint occupation by f a and d electrons at 
a site, the metallic state reappears. We have also cal- 
culated the corresponding density of states for different 
values of U and U / showing the metallic and insulating 
states as discussed above (Fig. 5 for 7i/i — nf 2 — 1/3). In 
Fig. 6 we show the corresponding phase diagram in U-Uf 
plane for nfi — u/2 = 1/3. 

The presence of the local U(l) gauge symmetry for 
each of the two localized levels have additional ad- 
vantages. As the conservation of f} a fi a renders the 
d— electron part diagonalisable (in a 'local' potential 
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FIG. 5: Density of states for different value of U and Uf at 
n /i = n Si = I/ 3 - 



of the /—electrons), the problem then becomes ex- 
actly solvable in the infinite dimensional limit where the 
d— electron self-energy is locate. Such a solution has 
been obtained by Brandt and Mielschii for the original 
Falicov-Kimball model. A straightforward generalization 
is possible in the present model as well. The local Green's 
function, in the limit of infinite spatial dimension, turns 
out to be 

G(u n ) =< (1 - ri/i)(l - n f2 ) > G (w„) 
+ < n/i(l -n/ 2 )+n /2 (l- n n ) > (Gq 1 (u n ) - U)~ x 

+ < 7i/in/ 2 > (Gq \uj n ) - 2Uy\(3) 

where Go(u) n ) = (iu> n + /i + <5(w„)) _1 with cu n — (2n + 
1)ttKT the Fermionic Matsubara frequency, 6(u) n ) the 
self-consistent, time-dependent generalized potential. In- 
serting the Dyson equation G ~ 1 (o;„) = G _1 (w n ) + S(w„) 
in the above, a cubic equation for the self-energy S(w„) 
can be written. If one further takes the limit of Uf —> oo, 
then the last term on the right side in Eqn.(3) above 
vanishes and writing rif = n/i + nf 2l one gets back the 
Brandt-Mielsch solution. Note that the Green's function 
appears identical to that of the coherent potential ap- 
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FIG. 6: Phase diagram showing metallic and insulating re- 
gions for nf 1 = rif 2 — 1/3, in U — Uf plane. The x-axis 
((7 = line) is trivially metallic. 




FIG. 7: The mean-field excitonic order parameter Af in the 
V-U plane for n/i = n/2 = 1/3. 
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FIG. 8: (colour online) The mean-field excitonic order pa- 
rameter for n/i = n/2 = 1/3 and U = 1 in the homogeneous 
(lower curve) and ordered (upper curve) ground states. Inset 
shows the orbital densities at A, B and C sublattice respec- 
tively. 



proximatior>i£ (as also alloy analogy, Hubbard-III), but 
the difference here is that < n,f a > are determined by 
the Green's function G(w„) through rij = f(Ef a ) and 
E fa = E fa - KT J2 u Jog(l - C/GqK)) (here f(x) is 
the Fermi function [l + exp(x — n) / KT]^ 1 ). In this limit, 
the self-energy is easily obtained in the case where the 
filling is < n/i >=< nf 2 >=< rid >= The choice of 
Hartree self-energy £ = U < n/ > evidently minimizes 
the ground state energy. The corresponding configura- 
tion is an ordered three-sublattice arrangement of d, /i 
and /2 electrons (with d, /i, f 2 electrons occupying A, B, 
C sublattices), in the limit U/ — > oo where the two types 
of /—electrons must not occupy the same site. In fact 
the argument, following Czycholl 1 — , can be extended (for 
the ground state at least) to the case where a small hy- 
bridization term like Vd\fi a breaks the local invariance 
of $ ct electron number, though the global U(l) symme- 
try is still extant. We turn to this situation now. 

There have been several studies on the FKM on a 
square lattice over the years with a hybridization term 
that mixes localized and extended electrons. An im- 
portant issue in this context was raised by Portengen 
et al., 20 regarding the spontaneous symmetry breaking 
(SSB) in such a model. A term Vd\fi. a would in- 
duce non-zero excitonic (or ferroelectric 2n ) averages like 
< <4/i,a > as long as V is non zero. It was shown^S from 
a mean-field theory calculation on a square lattice at half- 
filling that this average tends to a finite value even when 

V — > 0, leading to a spontaneously broken symmetry in 
the ground stated. Although the soft local dynamic fluc- 
tuations between f a and d electrons in the limit V — > 
cannot be treated properly in the static mean-field, we 
nevertheless undertake a similar mean-field analysis on 
the triangular lattice following previous wor k 19 ' 20 . As- 
suming a homogeneous ground state and excitonic mean- 
fields, we find (Fig. 7) that when U is finite, averages of 
the type A/ iQ =< d\fi a > do not vanish in the limit 

V — > 0, indicating an SSB similar to the situation on a 
square lattice 2 ^. Such a non-zero average would, then, 
imply a homogeneously mixed valent ground state. 

However, the assumption of a homogeneous ground 
state itself is called into question for the ranges of pa- 
rameters being studied. It is important to note that there 
could possibly be other ordered states with lower energy 
at half-filling and one would therefore need to first as- 
certain the correct ground state. Czycholl 19 had shown 
that on the square lattice, this is indeed the case and the 
right ground state is an ordered two-sublattice charge 
density state. He also showed that such a ground state 
does not support SSB and a homogeneous mixed valence. 
It is not clear a priori that on a triangular lattice, similar 
considerations would apply. Clearly our mean-field anal- 
ysis above shows that on a homogeneous ground state, 
the results of Portengen et al. hold. We, therefore, look 
for the ground state in the mean-field theory and work 
out the possibility of an SSB on a triangular lattice. 

Following the arguments in the infinite dimensional 



6 



limit above, we expect an ordered three-sublattice or- 
bital density state to appear at half-filling when rid = 
rif a = 1/3 even for small finite V. We set up the mean- 
field calculation to look for an orbitally ordered state and 
allow for inhomogeneous local order parameters. This 
involves calculating all the local order parameters in the 
real space iteratively for self-consistency. We describe the 
result in Fig. 8 where an inhomogeneous mean-field solu- 
tion with a three-sublattice orbital structure emerges as 
the lowest energy state. The orbital density order pa- 
rameter is shown in the inset of Fig. 8. The A sublat- 
tice has negligible /-occupancy and a finite d-electron 
density, while B (C) sublattice has predominantly /i 
(/2) occupation and a small d-occupation. Clearly as V 
rises, such a state becomes less stable and at a critical V 
the orbital order would melt (although mean-field anal- 
ysis may not be quite valid in that range) due to strong 
quantum fluctuations. On reducing V , we observe that in 
the limit of V — > 0, the excitonic order parameter (note 
A/ jQ , =< d\fi a > is same for a = 1,2 in the degener- 
ate limit considered) vanishes leading to an absence of 
SSB in this case. For comparison, we also show the be- 
havior of this order parameter in a homogeneous ground 
state. In this case the solution supports a non-vanishing 



expectation value even at V — > 0. 

In conclusion, we have studied an extended Falikov- 
Kimball model with two localized states on a triangular 
lattice. Such a model is found to reproduce orbitally or- 
dered insulating ground states as reported (or predicted) 
in experimental systems like Gdla^, NaTi02 etc. Fur- 
thermore, the disorder to order transitions have also been 
seen to accompany a metal-insulator transition. Many of 
these orbitally ordered ground states would induce spon- 
taneous lattice distortions due to the breaking of inver- 
sion symmetry and lead to ferroelectricity driven not by 
the conventional mechanisms, but by electronic correla- 
tions instead. We also investigate the possibility of a 
spontaneously broken symmetry in the form of an ex- 
citonic order paramete in this Hamiltonian and resolve 
that there is no SSB in the ground state on a triangular 
lattice. 
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